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The time process of transport on randomly evolving trees is investigated. By introducing the 
, notions of living and dead nodes, a model of random tree evolution is constructed which describes 

• the spreading in time of objects corresponding to nodes. It is assumed that ai t — the tree 

' consists of a single living node (root), from which the evolution may begin. At a certain time 

instant r > 0, the root produces u > living nodes connected by lines to the root which becomes 
dead at the moment of the offspring production. In the evolution process each of the new living 

. nodes evolves further like a root independently of the others. By using the methods of the age- 

, dependent branching processes we derive the joint distribution function of the numbers of living 

and dead nodes, and determine the correlation between these node numbers as a function of time. 
It is proved that the correlation function converges to VS/2 independently of the distributions of 1/ 
and T when gi — > 1 and t —> 00. Also analyzed are the stochastic properties of the end-nodes; and 
. the correlation between the numbers of living and dead end-nodes is shown to change its character 

O ' suddenly at the very beginning of the evolution process. The survival probability of random trees 

^ I is investigated and expressions are derived for this probability. 
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I. INTRODUCTION 



Recently there has been a growing interest in randomly evolving graphs since a large number of practically important 
problems can be analyzed in this way. In the last few years so many papers have been presented in this field that any 
list of references would be far from complete by any standards. It is fortunate that two outstanding review papers 
0, have been published in 2002. Thus the author of the present paper does not feel obliged to cite the very large 
^ ' amount of important references. In addition, one must mention the excellent book by Dorogovtsev and Mendes 
which gives a very interesting summary of the actual problems of the growth and structure of random networks. 

It has been long known that the random evolution of a tree ^ corresponds to a Galton- Watson process 4] which 
describes the formation of a population which starts ait — with one individual and in which at discrete time moments 
: t — 1,2, . . . each individual produces a random number of offspring with the same distribution independently of the 
Qv^ ' others. There are many interesting papers on this subject, statements of which are analyzed and generalized in Ref. 
T— I ! [3 . The purpose of this paper is to study the random tree evolution with continuous time parameter with the hope 
' that this may help to a better understanding of the nature of the age-dependent branc hing processes which play a 
^""^ _ significant role in many fields of physical and biological sciences 0,13, ^> IT^ H^ . The general theory of 
JiO ■ age-dependent branching processes is excellently summarized in Ref. 3 and in Ref^ therefore we decided to 

' focus only on the correlation and survival properties of the randomly evolving trees. 

I The paper is organized as follows. In Sec. II the exact description of the random process is given, while in 
^ ' Sec. Ill the derivations of the basic equations for the generating functions can be found. The properties of average 
characteristics (expectation values, variances, and correlation functions) are discussed in Sec. IV. Survival probability 

p. 

^ more detailed analysis of the problem see [iJI- 
O 



'"^ of trees is determined in Sec. V, and finally, the conclusions are summarized in Sec. VI. It is to mention that for 

o 



II. DESCRIPTION OF THE RANDOM PROCESS 



We will study the random evolution of trees which consist of living and dead nodes as well as lines connecting the 
' nodes. A living node is unstable, and after a certain time r may become dead. The lifetime r is a random variable, 
and T'jr e {t,t + dt)} = dT{t) denotes the probability of finding r in the time interval {t,t + dt). A dead node is 
inactive, it cannot change its state, and is unable to interact with the other nodes. 

The evolution process can be described as follows. Let us suppose that at time instance t — the tree consists of 
a single living node called the root and let iSo denote this state of the tree. At the end of its life the root creates 



A randomly evolving tree is a connected graph containing no cycles, and growing from a single node (root) according to well-defined 
rules. 
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V = 0, 1, 2, . . . new living nodes and immediately after that it becomes dead. The number v of the offspring is a 
random variable, and the probabilities V\v — — //j, fc = 0, 1,2, . . . are assumed to be known. The new living 
nodes are promptly connected to the dead one and each of them can evolve further like a root independently of the 
others. This means that the tree from its initial state may evolve further by two mutually exclusive ways, in which 
the root either 

• survives, and thus the tree remains in its initial state during the whole time interval (0,t), or 

• the root does die in an infinitcsimally small subinterval [t' ,t' + di') G (0,i) and produces j > new living 
nodes. Each of these nodes will evolve further in the remaining time interval (t', t) independently of the others 
like the root. 




FIG. 1: First step of the tree evolution. The new state arises from the initial state which consists of a single root. Four of the possible 
new states are illustrated in the figure. Living nodes are denoted by white circles and the dead ones by black circles, /o, /i, /2, /a, ■ • ■ are 
the probabilities of producing 0, 1, 2, 3, . . . nodes by one dying root. 

Figure. shows the initial state of a tree consisting of a single root and illustrates four of the possible states 
which may be produced by the dying root. The probabilities of producing 0, 1,2,3,... living nodes are denoted by 
/o, /l> /2, /3> ■ • ■ • 

Figure. |21 illustrates one of the realizations of a random tree at the moment i > 0. In Fig. [21 there are three white 
circles denoting living nodes which are capable of further evolution, and a large number of black circles denoting dead 
nodes, which are unable to produce new living nodes. It is to mention that each offspring is connected by a line to 
its dead precursor. If a dead node has fc outgoing lines it is called node of /cth out-degree. 

In order to simplify our consideration we will assume that the distribution function of the lifetime r of living nodes 
is exponential, i.e., 

dT[t) = e"'^* Q dt, (1) 

where l/Q is the expectation value of the lifetime. In this case the evolution becomes a Markovian process. It is 
assumed that the probabilities 

V{iy = k}^fk, A: = 0,1,2,... (2) 
are the same for all living nodes. By introducing the generating function 

oo 

= (3) 

fe=0 

we can define the factorial moments of v in the following way: 

r = l,2,... . 

z=l 



(Fqiz) 
d77 
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FIG. 2: One of the realizations of a random tree at tlic time moment f > 0. The nodes denoted by white circles are capable of further 
evolution, while those denoted by black cannot evolve and further the horizontal dashed lines indicate the time instances at which nodes 
were produced. 



For the expectation value and the variance of v we have 

E{z/} = gi and D^ji/} = 52 + 9i - 9? • 
For the sake of later use we cite the relation 

E{(z.-l)2} = D2M + (l-gi)2. 



(4) 



(5) 



In many cases we do not have to know the exact form of probabilities //j, A; € .E, it is enough to know only the first 
and second factorial moments of v. However, if we wish to obtain results which are, in some sense, exact then we will 
assume that fk = 0, Vfc > 2, i.e. we will use a quadratic generating function 



= fo + fiz + f2Z^ = 1 + - 1) + ^q2{z - if. 



(6) 



It seems to be useful to cite the following trivial relations: 



/o = 1 - gi + ^ ^2, fi=qi- 92, h = \ 12, 

which follow from the equations /o + /i + /2 = 1 , /i + 2/2 = ?! , and 2 /2 = 92 • Since /o , /i , and /2 are non- negative 
real numbers smaller than 1, and their sum is equal to 1, the possible values of g'l and q2 are restricted. 
As a source of simplification, sometimes we use the geometric distribution 



= k} 



1 



91 



1 + qi VI + 91, 

which is defined by the single parameter qi only. The generating function is expressed by 

1 



'^^^'^ = l + (l-.)9i^ 



(7) 



and we find immediately that E{j/} = qi and D^{j/} = qi(l + gi), i.e., (72 = 2gJ. 

In order to characterize the random process of tree evolution we introduce the random functions ^ie{t) and Hdit) 
which are the momentary numbers of the living and dead nodes, respectively, at time instant t > 0. It seems to 
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be worthwhile to define the random function fi{t) = fJ.e{t) + fJ'd{t) which gives the total number of nodes, i.e., the 
momentary value of the tree size at t > 0. 

The nodes can be sorted according to the number of outgoing lines. Let us denote by xei't, k) the number of living, 
while by Xdit,k) the number of dead nodes with k outgoing lines, at the time moment t. A node not having any 
outgoing lines is called an end-node. Obviously, all living nodes are end-nodes, i.e., Xi{tjO) — ^ii{t) and xeit,k) = 0, 
if fc > 0. Clearly, 

oo 

fJ-d{t) ^^Xd{t,k). 

k=0 

To simplify the notation, we introduce the following random functions: 

Xeit) = Xiit,Q) and Xd{t) = Xd{t,0). 

The total number of end- nodes is given by x{t) = X^(0 ~^Xd{t)- In the next section we shall derive the basic equations 
for the generating functions of probabilities of the random events — ng, ^ld{t) = 'T-d}, {/^(O — and 

{xt{i) — T^i^ Xd{i) = nd}^ {x{t) — provided that at t = the tree was in the state Sq. 

III. DERIVATION OF THE BASIC EQUATIONS 

In this section we would like to deal with equations for the generating functions 

5(^''^)(t,z,,z,)=E{zf(*) and z) = E {z^^ } , (8) 

as well as for 

#'^)(i,z,,z,)=E{zf(*)zf(*)} and z) = E {z>^W } . (9) 

It is seen immediately that the generating functions g^^\t, z) = E {z^**^*^} and g^'^\t, z) = E |z^<i(*)} can be obtained 
from © since 

g(''^\t, z, = z, Zrf = 1) = z), g<^''''\t, z, = 1, z, = z) = g^^\t, z), 

and similarly, gi^\t, z) and gi'^\t, z) from lO, so it is clear that we need equations only for the generating functions 

g'-'^-'^\t,ze,Zd), g{t,z) and gl^''^\t, ztZd), ge{t,z). 

It is evident that these equations belong to the class of equations corresponding to the age-dependent branching 
processes. 

A. Joint distribution of the numbers of living and dead nodes 

We shall define now the probability 

V{Mt) = ^H,^id{t) = nd\So} = p^^^'^\t,ni,nd), (10) 

where ni and Ud are non-negative integers. It is clear that ng, Ud) gives the probability that at the time instant 

t > the number of living and that of dead nodes in a randomly evolving tree are equal to and n^, respectively, 
provided that at t = the tree was in the state 5o. By exploiting the rules of random evolution described in Sec. II 
and applying the Theorem 7.1 of Ref. (3| it can be shown that 

p^^^'^\t,ni,nd) = e~^* 5„^,o+ 
Q ^ e-«(*-*') |/o ^„„o V'^''''Ht' ,nt,nd\k)^ dt, 



where 



i-nt,k=ni nd,iA hnd.i—nd - 1 j = l 



Introducing the probabihty generating function 



oo oo 



n^— Ufi—O 

and taking into account the generating function defined by Eq. ((SJ, we obtain the integral equation in the form 



Qt 



dt\ 



which according to the Theorem 11.1 of Rcf. 4] can be rewritten as a nonhnear differential equation 



dt 



with initial condition 



-Q g^'^^\t, ze, zd) + ZdQq g^^'^^'it, zt, z^) 



lim g^^''^\t,Zi,Zd) = z^. 



(12) 



(13) 



1 . Remarks 



By using the relations g^^''^\t, zi — z,Zd = 1) ~ g^^\t, z), g'-^''^\t, zi — l,Zd ~ z) ^ g^'^Ht, z) from \\2\ we obtain 
the following two equations: 



g'^'\t,z)=e-'^' z + Q 



-Qt 



{t',z) 



dt' 



(14) 



-,{d)( 



(t, z) = e-'?* + zQ e-«(*-*') q U''^ (f , z) 



dt'. 



It is an elementary task to prove that the integral equation H14|) is equivalent to the differential equation 



dg'^'Kt^z) 
dt 



= -Qg^'Ht,z) + Qq g^'>it,z) 



with the initial condition limt|o g^^^tT^) = z, while Eg. ljlSII to 



dt 



Qg(^^{t,z) + zQq g(''\t,z) 



with the initial condition limtjo g'^'^H^T^) — ^■ 

Finally, it can be easily shown that g{t, z) — E{z'^(*)} satisfies the integral equation 

g{t,z)^e-^' z + zQ f e-^^*'*'^ q[g{t' , z)] dt' . 
Jo 

For the sake of completeness we derive from Eq. (|18|l the equivalent differential equation 

dg{t,z) 



dt 



~Q g{t,z),+Qz q[g{t, z)] 



the initial condition of which is limtjo g{t^ z) = z. 

It can be shown jl5| that each of the equations for the generating functions has just one regular solution. 



(15) 



(16) 



(17) 



(18) 
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B. Joint distribution of the numbers of living and dead end- nodes 

Now we would like to determine the joint distribution of xe{t) Xd(t), i.e., the probability 

V{xi{t)=ne, Xd{t) ^ nd\So} ^ pi''''Ht,ne,nd). 

It is evident that 



+Q / e 

Jo 



-Q{t-t') 



k=l 



dt', 



where 



"*,lH [-n«,fc=nf nd,iH Vnd,k=nd j = l 

Simple calculations show that the generating function 

oo oo 

satisfies the equation 



m—Q Ud—O 



(19) 



which will be used for the determination of the correlation between the random variables xe{t) Xd(t). For the 
sake of brevity the corresponding differential equation is omitted. 



1. Remark 



There is no need to write down the equation for generating function of x/i't)j since gi\t,z) = g^{t) which is 
the solution of Eq. H14() . The equation for the generating function gi'''\t,z), however, one can obtain from l|19() by 
substitutions zi = I and Zd = z, and so one has 



gi'Ht, z) = e-^' - (1 - z)h{l - e-<^^) + Q 



-Qt\ 



dt'. 



(20) 



In order to gain an insight into the interplay between the living and dead end-nodes it seems to be useful to calculate 
the probability distribution of the random function x{t) =^ Xei^) + Xd{t)- The equation for the generating function is 
easily obtained as 



,{t,z) = J2'^{xit) = n\So} z" = 5^Pe(i,n) 



where Pe(t,n) is the probability that the number of all end-nodes at t > is equal to n provided that at i = the 
tree was in the state Sq. Applying the usual argumentations one has 



ge(t,^) = e-«*z-(l-z)/o(l-e-'?*) + g / e-Q('-''k[g„it' , z)] dt' . 

Jo 



(21) 



IV. AVERAGE CHARACTERISTICS 



A. Correlations between the numbers of living and dead nodes 



The study of how the random functions iJ,e{t) and fid{t) are related involves the analysis of correlation between 
them. First we should calculate the covariance function 



Cov{neit)nd{t)} = B{iieit)nd{t)} - E{ni{t)}E{^d{t)}, 
and then the correlation function 



r{i,d)U n r,\- Co-v{^t,{t)^ld{t)} 

'Kt,qi,q2) — — — , , — T-TY- 
D{iie{t)}Ti{iid{t)} 



(22) 



(23) 



As is apparent, we must determine the expectations and variances of iJ,e{t) and iJ.d{t). Taking into account the 
relationships 



and 



E{Md(<)} 



dz 



\ dz 



2=1 



mf(t) 



after simple algebra we find 

mf^(t)=e-"* and mf\t) 



^ [l-e-«*], if^i^l, 

Qt, if qi = 1, 



(24) 



where a = {1 — qi)Q. It is here to mention that the tree evolution is called subcritical, if gi < 1, critical, if gi = 1, 
and supercritical, if gi > 1. 
The formulas 



D^{/x^(i)} = (i) + w>^'{t) 1 - m!(>{t) 



and 



show that we must determine the second factorial moments 



(25) 



(26) 



After elementary calculations we find that 



and 7712"^^ {t) 



2=1 



it) 



_^ [1 _ e-"*] e-«*, ifgi^l, 



and 



{Qt? + \ 92 {Qtf, if qi = 1, 



92 Qt, 



if 91 = 1, 



(27) 



m('^'(t) 



(28) 



where 



and 



(29) 



id) _ 92 

By substituting the first and second factorial moments of Eq. H25|) for formulas (|24|) and H27|l , we obtain that 

T>^MQt, ifgi = l, 

where E{(zy — 1)^} = (1— (71)^ + D^{i/}.^ Applying similar substitution in H26() . if (71 7^ 1, then we arrive at 

n^fiaim = vi"^ [1 - e-"*] + v!f^ Qt e^"* + «f e""* [1 - e^"*], (30) 

where 

Jd) _ DH^} „(<i) 



If (71 = 1, i.e., the evolution is critical, then 



D^z.} (rf) _ E{(^-l)n 



(1-91) 



1 



2' ^3 



(1-91)^ 



(31) 



One can see that the limit values of the variance of the number of dead nodes are given by the formula 

. (1-91)^' 

hm D^{/Xd(i)} 

t — ^oo 



14 , if qi < 1, 



00, if 51 > 1. 

In order to calculate the covariance function wc should solve also the equation determining the mixed moment 



E{fie{t)^ld{t)} = 
The equation can be written in the form 



dzidzd 



(32) 



— Zd — 1 



By using Eq. H24() for m,f\t) and m'"^\t), and introducing the notation a = (1 — gi) Q, it can be proved that 



Q I e 
/o 



91 ml'\i') + 92mf (f) m[''\t') 



dt'. 



(33) 



^ e""* - 52 e-"* (1 - e^"*), if gi ^ 1, 



Qi+iD2{j.}(Qi)2, if 91 = 1 

is the unique solution of the integral equation H31|l. After a brief calculation we have 

e-"* (1-e-"*), ifqi ^1, 



Cov{^^(i)^d(i)} = 



1-91 



1 -L D i'^}. 



(34) 



if qi = 1. 



^ In the case of subcritical evolution the variance has a maximum at tmax = log2/a which, as seen, is independent of the form of the 
distribution of p. 
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time in Q"^ units 

FIG. 3: Three curves illustrating the time dependence of the covariance Cov{iJ,£{t)fi^(t)} in the case of subcritical random evolution 
with the assumption that u is of geometric distribution. The curves a, b, and c correspond to the values of qi = 0.96, 0.97 and 0.98, 
respectively. 



The dependence of the covariance function Cov{^g{t) fid{t)} on the time parameter Qt at three different values of 
qi has been calculated by assuming that v is of geometric distribution. The three curves a, b, and c corresponding to 
the values qi = 0.96, 0.97, and 0.98, respectively, are plotted in Fig. |31 One may observe that each of the covariance 
curves vs time has a well defined maximum. The time parameter which belongs to the maximum of the covariance 
function depends rather sensitively on qi. 

Now, we would like to return to the study of the correlation function defined by (|23|l . It is evident that the numbers 
of living and dead nodes are correlated, but it is worthwhile to note that the time dependence of the correlation 
function C'-^'''-' (t, gi, 92) has some specific features, which follow from properties of the covariance function and the 
variances of ne{t) and fidit). For example, one can show that in the case of subcritical evolution the correlation 
function, after a sharp increase, converges to zero when Qt —^ cxd, while in the case of supercritical evolution it rapidly 
reaches the value 1. This behavior is illustrated in Fig. 2|for the case of ly following a geometric distribution. The 
curves a, b, and c belong to values 0.95, 1, and 1.05 of qi, respectively. It seems to be interesting to present some 
curves reflecting the dependence of C'^^''*^ (t, gi, 92) on qi at different time instants. In Fig. |31we can see that the 
transition of the correlation function from to 1 becomes sharper and sharper as Qt increases. The curves a and b 
belong to the time parameters Qt = 100 and Qt = 500, respectively, while curve c shows how the correlation function 
depends on qi at Qt = 00. 

The most surprising result can be obtained when the random evolution is critical, i.e., qi = 1. Then we have 



1 



which shows that the limit value of the correlation 

hm C<''''\t,l,q,) = 

Qt — *oo 



B^iy} [Qtf 
V3 



-1/2 



(35) 



(36) 



is independent of the distributions of v and t. This very important limit law expresses the fact that the relation 
between the numbers of living and dead nodes near the critical state, and at time moments sufficiently far from the 
beginning of the process, becomes almost independent of rules controlling the evolution. This independence can be 
attributed to the appearance of the universal properties of the evolution. 
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time in Q"^ units 



100 



FIG. 4: Time dependence of the correlation between the numbers living and dead nodes if u is of geometric distribution. The curves a, 
b, and c belong to gi values 0.95, 1, and 1.05, respectively. 
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FIG. 5: Dependence of the correlation between the numbers of living and dead nodes on qi at time instants Qt = 100, (a), Qt = 500, (b) 
and Qt = oo, (c). 



B. Correlations between the numbers of living and dead end-nodes 

The dynamics of the random tree evolution is determined by the end-nodes. It is evident that living end-nodes 
are responsible for the development of a tree, while dead end-nodes represent points where the tree is not capable of 
any variation. In order to have an insight into the dynamics of the tree evolution it is useful to study the correlation 
function of xeit) and Xd{t)- First we calculate the covariance function 



Cov{xe{t)xdm = E{xe{t)xdm - E{x^(t)}E{xd(i)}, 



(37) 
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and then the correlation function 



D{xHO}D{xd(i)} 



(38) 



Since — Mf(^)7 we have to deal only with the determination of the variance D^jxd}- Omitting the details of 

calculations, if qi 7^ 1, we have 



(1-91)^ 



g2 

1-91 



1-91 



-1) (1 - e-"*) - 2g2Qi' 



92 



,1-91 

while, if gi = 1, we obtain the following simple formula: 



1 (1 



D'{Xd(t)} = /oQM l + :792(9i) 



(39) 



(40) 



It is worthwhile to remark that if the tree evolution is subcritical, then and only then does the variance of the number 
of dead end-nodes converge to a finite limit value given by 



lim V>\xm 



(l-(7i)2 Vl-91 



92 



1-91 
/o 



- 1 



As regards the covariance function from the generating function gi^''^\t, Zi, Zd) defined by (|19|l we should derive the 
second mixed moment 



m['f{t\e) = 



d^gi^''^\t,Zi,Zd) 



dzt dzd 



which is determined by the following integral equation: 

nt 



m[^f\t\e) = qi 

Since 

we obtain the solution in the form: 







Q(t-t') m['f{t'\e) dt' + q2Q / e-'?^*-*') m['\t'\e) m[''\t'\e) dt' . 



(41) 



and 



1-91 ' 



/o 



Co-v{xi{t)xd(t)} = 



1-91 



1(1 



if 917^1, 



if = 1. 



(42) 



As is to be seen in Fig. El the covariance of the numbers of living and dead end-nodes changes its sign at the very 
beginning of the evolution. The negative covariance just after the start is a consequence of the simple fact that the 
root of the tree can be but a living end- node at the time instant t = Q. The calculations refer to the case when v is 
of geometric distribution; this however, does not affect the generality of the conclusion. 

For the sake of completeness we calculated the time dependence of the correlation function of the numbers of living 
and dead end-nodes at values of qi — 0.95, 1, 1.05 assuming again that v is of geometric distribution. As expected. 
Fig. 13 shows convincingly the change realizing in the character of the correlation near the beginning of the evolution. 
There is no need to write down the lengthy expression of the correlation function but it seems to be useful to present 
it when the tree evolution is critical. In this case one obtains 



2 Y 92 V3 + 92(Qt)2 



(43) 



It is interesting to note that the correlation in critical evolution reaches its asymptotic value ^3/o/2 very soon. (If 
/o = 0.5, then cf''''\oo) ~ 0.612 and this is what is seen in Fig. [7|) 
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FIG. 6: Dependence of the covciriance of the numbers of hving and dead end- nodes on Qt at different values of qi. 
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FIG. 7: Dependence of the correlation between the numbers of living and dead end-nodes on Qt at different values of qi . 

V. SURVIVAL PROBABILITY 

A. General considerations 

It is obvious that the evolution of a random tree will stop at the time instant 9 which satisfies the equation i^i{0) = 
with probability 1. The random variable 6 is called the lifetime of the tree. In order to determine its distribution 
function 

V{e < t\So} = m, (44) 

one must recognize that the probability V{f-ii{f) = 0\Sq} = p*-^^ {t, 0) to find zero living node at the time moment t > 
in a tree is the same as the probability that the lifetime 9 of that tree is not larger than t > 0, therefore, one can 
write V{9 < t\So} = V{ni{t) = 0|5o}, i.e., 

L{t) = p^^\t, 0) = lim g^^\t, z). (45) 

It is clear that if < < ^2 then L{ti) < L{t2), i.e., L{t) is a nondecreasing function of its argument, hence the limit 
relation 

max L{t) = lim L{t) = < 1 (46) 

0<t<cso t—*oo 
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must hold. We will call the quantity Loo extinction probability, and will now prove the following statement: if qi < 1, 
i.e., the random evolution is not supercritical, then L^o = 1, while if gi > 1, i.e., the evolution is supercritical, then 
Loo is equal to the non-negative, single root of the function ■0(2/) = q{y) — y, y & [0, 1]. (It is evident that 0(1) = 0.) 

For the proof we exploit the fundamental property of the generating function g^^') {t, z) which is expressed by the 
equation g'^^\t + tt, z) = g^^^ \t, g'^^^u, z)] . Applying the relation (|46|l we have L{t + u) = g^^-* [t, L{u)] and since 

lim L{t + u) = lim L{u) ~ Loa, 



we can write for every t > that L^c — g^^H'^j Lao)- Plugging g'^^(i, L^o) into the equation 



dg^'\t,. 



g(^)(t,z) -Qg^'>it,z), 



we obtain 



g (Loo) 



0. 



(47) 



(48) 



Considering that q(y) is a probability generating function, i.e., limj,|i q{y) — 1, then according to a well-known 
theorem of generating functions, it is clear that ii qi > 1 then Eq. (|48|l besides the trivial fixed-point 1 must have 
also another, smaller than 1, non- negative fixed-point L^o, and this is what we wanted to prove. 
Now, we would like to define the survival probability. Obviously, 



S{x) = 1 - L{x), 



(49) 



is the probability to find the tree at the time instant x = Qt in a living state, and therefore, S{x) can be called the 
survival probability. By taking into account the properties of L{x) one obtains 



0, if qi < 1, 

5oo = 1 - Loo, if qi > 1- 
Clearly, from Eqs. (|47|l . (|45|l . and (|49|l one can derive the equation 



lim S{x) 

Qt^oo 



^ = -g(l-5)-5 + l. 



which has the solution 



at the initial condition 5(0) — 1. 



x{S) 



dy 



S{x 



) g(l - J/) +y- 1' 



(50) 



(51) 



(52) 



B. Calculations in the case of a quadratic q{z) function 

In order to demonstrate the characteristic features of the lifetime of random trees we will determine the survival 
probability in the case when the generating function of the offspring number is a quadratic expression defined by 
©. By using Eq. (|FT|l we obtain 



- = -il-q,)S--q,S , 
and taking into account the initial condition 5(0) = 1 we have the solution in the form 



(53) 



Six) 



92 



2(l-9l) 



(l-e-(i-«i)^) 



2+92 X ' 



2 ^ 

92 



l + (l_2 2i^)e-(«i-i)=^ 

^ 92 ' 



, if gi < 1, 
if qi = 1, 

1 

if qi > 1. 



(54) 



14 




FIG. 8: Dependence of the probability density function of the lifetime on x = Qt. 



The density function of the hfetime measured in units of Q ^ can be calculated by using the relation 

dL{x) dS{x) 



l{x) 



dx 



dx 



It is elementary to show that the density function is decreasing monotonously from ^(0) — fo to zero. In Fig. (Hjone 
can see the density function curves versus time x — Qt for subcritical, critical, and supercritical trees. 
From H54|l one can see that 



lim S (x) 

x—Qt — >oo 



0, if<Zi<l, 
5^ = 2 21^, ifgi>l. 



(55) 



For more detailed analysis of asymptotic values for S{x) see Ref. [l^ . 

Finally, in this section let us calculate the expectation and the variance of the tree lifetime 9 in units of . For 
this purpose it seems to be useful to determine the characteristic function of the random variable Q9. Since the 
moments E{((36')-'}, j = 1, 2, . . . do not exist if qi > 1, the calculations are restricted to the case when qi < 1. One 
can write 



POO pi 

^{iu) = E{e-"^^} = / e""" dL{x) = / , 
Jo Jo 



-ujx{y) 



where is a complex number with Reuj > 0. Performing the substitution 

1 



l~qi 



l + 2i^ 

In y 



y + 2 



92 



one obtains the characteristic function 



where 



/? = (1 - <7i)-i and 7 



/" 


y 


Jo 


.y + 1. 






7 = 






92 



w/3 



dy, 



(56) 



qi < 1. 



From the characteristic functions ip{uj) both the expectation and the variance of the lifetime can be easily calculated 
when the tree evolution is subcritical. For the expectation value one obtains 



E{Q0} = - 



/ d(p{u}) 
\ du! 



= — In I 1 

92 



1 92 



2 1-91 



(57) 



li qi ^ 1, then the expectation value diverges as In (1 — gi) ^. 



15 



For the calculation of the variance we need the second moment of which can be immediately obtained from the 
characteristic function. Omitting the details of the calculation we may write 



-) 



1 92 



2 1 



1\ 



In 1 



1 

2 r 



92 



9i 



n 2 



92 (1 - 9i) 



Li. - 



92 



2 1 



91 



(58) 



where Li2{u) = ^'^/^^ ^he so called Jonquicre's function. When qi is approaching 1 from below, the 

fluctuation of the tree lifetime becomes unlimitedly large, and so in the vicinity of the critical evolution the average 
lifetime loses almost completely its information content. 



VI. CONCLUSIONS 



Let us summarize the main results of the paper. By introducing the notions of living and dead nodes, a model of 
random tree evolution with continuous time parameter has been constructed. The model describes the spreading in 
time of abstract objects which correspond to nodes. It is to mention that the process of the random tree evolution 
analyzed in this paper belongs to the family of the age-dependent branching processes. 

In order to characterize the evolution process, essentially two basic random functions iJ.e{t) and Hdit) have been 
used. The first one is the momentary number of living nodes, while the second one is that of dead nodes at the time 
instant t > 0. By the assumption that the evolution process controlled by the lifetime r and the offspring number u 
of living nodes, exact equations have been derived for the generating functions of the probabilities 

V{ne{t) = ne, iid{t) = nd\So], and V{^,{t) = iii{t) + = n\So}, 

where Sq indicates that at t = the tree consisted of a single living node only. To complete the description of the 

tree dynamics the notion of end- nodes has also been defined. It is remarkable that the average lifetime of living nodes 
has a role in scaling the time in the generating function equations, only. 

The time dependence of the expectation value of the number of living nodes has been analyzed, and as expected, 

there exist three completely different types of evolution depending on the average number qi of offspring produced 
by a single living node, li qi < 1, then the evolution is subcritical, if gi — 1, then it is critical, while in the case of 
g'l > 1 is supercritical. 

A specific property of the tree evolution has been discovered, namely, it has been proved that the correlation between 
the numbers of living and dead nodes decreases to zero in subcritical and increases to 1 in supercritical evolution, if 
the time parameter Qt tends to infinity, but in the case of exactly critical evolution it converges to a fixed value \/3/2 
which is free of the process parameters. 

The stochastic properties of the end-nodes have also been analyzed, and it has been shown that the correlation 
between the numbers of living and dead end-nodes changes its character rather suddenly at the very beginning of the 
evolution process. It has been proved that the correlation function is negative just after the start of the process, but 
it becomes positive very sharply after elapsing a relatively short time. This behavior is the consequence of the fact 
that the root of the tree can be but a living end- node at the time moment t = 0. 

For the sake of better understanding of the evolution the survival probability of random trees has been investigated, 
and exact expressions have been derived for this probability in the cases of subcritical, critical, and supercritical 
evolutions when the generating function of the offspring number u was quadratic. 
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